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Renormalisation Scale Dependencies in Dijet Production 

at HERA 
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Abstract: Different choices of the renormahsation scale (uren) can be used to de- 
scribe hard scattering processes with two jets at large transverse momentum in 
deep-inelastic scattering at HERA by fixed order perturbative QCD calculations. 
For leading (LO) and next-to-leading order (NLO) calculations the simplest choices, 



0\ • Q^, the virtuality of the incoming photon, and the mean squared transverse mo 



a^ : ienta of the tw; jets E^ are studied in different kinematic regimes. It is found 



that only if both Q^ and E^ are large, the NLO calculation is stable with respect 
^ I to numerical variations of firen while the LO calculation strongly depends on it. If 

only one of the two scales, either Q^ or E^, is large, the NLO is more stable than 
rg ' the LO calculation, but exhibits nevertheless a strong residual scale dependence. 

When both scales are relatively small, large scale dependencies are found in both 
cases. Moreover, large differences between the LO and NLO calculation are found. 
Generally, the use of E^ as renormalisation scale is favoured over Q^, since scale 
j^ I dependencies are less pronounced and NLO corrections are smaller. 

O 

0\ ■ 1 Introduction 

0\ 

'^ 

p i' HERA colliding 27.5 GeV positrons on 820 GeV proton offers an ideal testing ground for 

A ■ perturbative QCD (pQCD) in deep-inelastic scattering (DIS). The centre of mass energy of 

O . ^ys ~ 300 GeV leads to a large phase space for hadron production and to the possibility to 

' observe collimated sprays of hadrons - called jets - in the hadronic final state. They are the 

experimentally accessible signs of a hard scattering process and relate the fundamental objects 

/\ ' of pQCD, the unobservable quarks and gluons, to the measurable hadronic final state. 

5t , At HERA, events with two jets in the central part of the detector can be produced in quark 

ill ~^ 19) or gluon {gj — > qq) initiated hard subprocesses (see Fig. p. In fixed order pQCD 
the dijet cross section can be written as: 



d (Tdijet Y^ „/ 2 N [^ dC, 



^^^Q2 E<(/^ren) /^ J Cn{x/ ^, fi^^, f^%,, ••.) ■ PDF{^, f,}j (l) 

The variables x and Q^ are the usual variables to inclusively describe DIS (see Fig. |I| for defini- 
tion). The (non-perturbative) universal parton density functions of the incoming proton PDF 
are evaluated at the factorisation scale fi'jac- They further depend on ,^ = a; (1 + s/Q"^), where 
V§ is the invariant mass of the hard partonic system. The variable ^ can be interpreted as the 
fractional momentum of the parton initiating the hard subprocess with respect to the proton 
momentum. In leading order, where s is the invariant mass of the dijet system, ^ can be di- 
rectly measured once the dijet system is tagged. The coefficient functions C„ can be calculated 
in pQCD as power series expansions in the strong coupling constant as{^ren) probed at the 
squared renormalisation scale fi^en- 

^ Contribution to the Proceedings of the DESY Workshop 1998/99 on Monte Carlo Generators for HERA 
Physics. 





DIS variables: 

x = QV(2P-7), y = Qy{sx) 
W'^ = Q^ {l-x)/x 

Xg/p = Xq/p = ^ = X (Q"^ + S)/Q'^ 



e^p - 



Figure 1: Feynman diagrams for the production of dijet events to first order of Ug in 
collisions. 7 (P) denotes the four-momentum of the photon (proton), ji and J2 are the four- 
momenta of the jets associated to the hard subprocess. 

The dijet cross section including higher order parton emissions are difficult to calculate using 
exact expressions for the coefficient functions. Since the phase space integrals cannot be solved 
analytically, numerical methods have to be applied. Several Monte Carlo integration programs 
are available to compute jet cross sections in next-to- leading order (NLO) of ^^(MEPJET ||I], ^, 
JETVIP [|, DISENT [g § and DISASTER g). They have been carefully compared in this 
workshop [0. They allow any jet definition scheme and arbitrary experimental cuts to be 
analysed. In this article all results are based on the DISENT program. NLO (LO) calculations 
are performed using the CTEQ4M (CTEQ4L) parton density functions [If. The value of as 
is taken from the parton density functions. In CTEQ4L as{Mz) = 0.131 and in CTEQ4M 
as{Mz) = 0.116 is used. The factorisation scale has been set to fi'j^^ = Q^. The dependencies 
of the results presented here on the factorisation scales are small and can be neglected. 



2 Choice of the Renormahsation Scale in DIS 



Dijet production in 
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Figure 2: Diagram for the production of dijet 
events in the Breit or hems frame in DIS. 



In pQCD calculations to all orders of a^ the 
dijet cross section does not depend on the un- 
physical quantities firen and fifac- However, in 
a calculation to finite order a residual scale de- 
pendence is expected. The size of the scale 
dependence can be used to estimate the higher 
\ proton order contributions which have been neglected 

in the calculation. It is expected that they 

are smallest when the perturbative expansion 
is performed at the "natural" scale. In DIS 
jet production is a two-scale problem where 
it is not a priori clear at which scale a.^ and 
the parton density function should be probed. 
In a frame0, where the virtual photon and the 
proton collide head-on, it is obvious that two 
scales can characterise the hardness of the pro- 



/ 



cess: Q , the virtuality of the incoming photon, as used in the QCD analysis of the inclusive 



^e.g. the hadronic centre of mass (hems) frame defined by 7 + P = or the Breit frame defined by 
7 + 2a; P = 0. In the laboratory frame and in the Breit frame the proton moves towards the +z direction. 



DIS cross section, or the mean squared transverse momenta of the two jets E^ as used in 
hadron-hadron coUisions. In this article the behaviour of the dijet cross section in different 
kinematic regimes for both Q^ and E^ as nominal renormalisation scales is investigated. 

To decide whether a physical process is probed at its natural scale and whether the pQCD 
prediction can be expected to be stable, the following criteria could give some guidance. The 
NLO corrections, i.e. the ratio of the NLO to the LO cross section, the so-called K-factor, 
should not be too large, since in this case it is conceivable that corrections of the next-to-next- 
to-leading order can be large. The residual scale dependence, i.e. the variation of the cross 
section when changing the nominal squared renormalisation scale {fJ'nom) by about an order of 



magnitude using a scale factor 1/4 < ^ren < 4, should be small (/x, 
the NLO calculation should depend less on variations of ^ren- 



2 
ren 



^ren- fJ^ioJ- Morcover, 



3 Dijet Rates and Cut Scenario 
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Figure 3: Dijet rate R2 = N2/NDIS as a function of the threshold difference A of the first and 
second jet. Jets are defined by the cone algorithm (Rcone = ^) in the hadronic centre of mass 
system with Et2 > 5 GeV and En > 5 + A GeV. Shown are HI data (points) for 11 < (^^ < 15 
GeV^ (a,h) and 50 < Q^ < 100 GeV^ (c,d) together with a NLO QCD calculation (line) with 
Q^ (a,c) and the mean E^ (b,d) as renormalisation scale fJ^nom- '^^e shaded band indicates the 
variation of the NLO prediction for 4 • finom o^nd 1/4 ■ n^om- 



Before studying the scale dependencies of the dijet cross section, it is first indispensable 
to find a cut scenario where pQCD can describe the data. For dijet production in DIS this 
has been a problem at HERA. For a long time agreement with NLO calculations could only 
be achieved in very restricted phase space regions (see e.g. [^j). Especially in the region of 
low Q^ < 50 GeV^, NLO calculations failed to describe the data. Meanwhile, it has been 
understood that this discrepancy occurs when both jets are required to be above the same Et 
threshold. As A defined as the required Et difference of the threshold of the two leading jets, 
approaches 0, a fixed order calculation gets infra-red sensitive |jlO|, 111], O]. In this region a 



reliable prediction of the jet cross section is not possible with a fixed order calculation. 

This has been demonstrated in an analysis by the HI collaboration ||T^, where R2 = 
N2/N£)js, the ratio of the two jet cross section to the inclusive DIS cross section, has been 



measured in the regioif] 5 < Q^ < iqo GeV^ 156° < 6^1 < 173", E^i > 11 GeV and y > 0.05. 



Jets are defined in the hadronic centre of mass system by the cone algorithm |1J] with a radius 
Rcone = 1- The Et of the second jet is required to be Et2 > 5 GeV and the Et of the first 
jet Eti > 5 + a GeV, where A is set to or 2 GeV. The difference of their pseudo-rapiditiesQ 
should be \ri*\ < 2. 
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Figure 4: Dijet rate R2 = N2/NDIS o,s a function of Q^ . Jets are defined by the cone algorithm 
(Rcone = ly' in thc hadronic centre of mass system with Et2 > 5 GeV and En > 7 GeV (a,c) 
or Et2 + Eti > 13 GeV (h,d). Shown are HI data (points) together with a NLO QCD calcu- 
lation (line) with Q^ (c,d) and the mean E^ (a,b) as renormalisation scale finom- "^^e shaded 
band indicates the variation of the NLO prediction for 4 ■ /i^o^ and 1/4 



r^nom ' 



The data points for 11 < Q^ < 15 GeV^ (a,b) and 50 < Q^ < iqo GeV^ (c,d) are shown 
in Fig. |. As expected, in the data the dijet rate is higher at A = GeV than it is at 
A = 2 GeV. In the NLO calculation the expected rise of the dijet rate towards decreasing 
A is also visible. However, around A k, 0.5 GeV the calculated dijet rate turns around and 
falls down until A ?^ GeV is reached. Around this point the NLO cross section falls with 
an infinite slope |T^, but remains finite. In this region the E^ of the jets are approximately 
equal and no phase space is available to emit a third real parton. This leads to an incomplete 
cancellation between real and virtual corrections at the threshold and makes a fixed order 
calculation unpredictive. A resummation of higher orders is necessary at this phase space 
point. Since the virtual corrections give a negative contribution, the dijet rate artificially drops 
down in the NLO calculation. The problem only occurs for events containing two jets with 
balanced momenta. In such a configuration the emission of a real third parton would lead to a 
configuration where one parton is above the Et threshold and the other below. If the thresholds 
of the two jets are different, the event has either balanced jet momenta or has jets with one 
above and the other below the thresholds. Since these conditions are never fulfilled at the same 
time, no problem occurs. 

In Fig. ^ one can see that the point A = 2 GeV is well described by the calculation while 



^ Of^i {Eel) is thc polar angle (energy) of the scattered electron. 

^ The pseudo-rapidity is defined as 77 = —In tan 0/2. In the laboratory and in the Breit frame the proton 



moves into the +z direction. 



for A = GeV the calculation is below the data. This behaviour is independent of the used 
renormalisation scale. However, when Q"^ is used a renormalisation scale a much larger scale 
dependence is found. In Fig. ^ and ^c it is demonstrated that the NLO calculation agrees well 
with the data for 5 < Q^ < 100 GeV^, if A = GeV is avoided. Hadronisation corrections 
would lower the NLO prediction by about 10 — 20%. Only in the case of the lowest Q^ bin the 
NLO calculation falls below the data. The difficult phase space region can not only be avoided 
by a cut on A, but also by a cut on the sum of the two jets, e.g. Eti + -Et2 > 13 GeV (see 
Fig. ^3 and ^d). For both cut scenarios the Q^ scale exhibits a much larger scale dependence 
than the E"^ scale (where E^ = {Exi + -Et2)/2). The variation of the NLO calculation when 
changing the renormalisation scale by a factor of 4 up and down is illustrated as grey band. 
The central value for a scale factor of C,ren = 1 is given as solid line. For both choices of scales, 
Q^ and E^, the cut on the sum of the two jets leads to a slightly smaller scale dependence than 
the A = 2 GeV cut. 

4 Dijet Cross Sections and their NLO Behaviour 
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To investigate the interplay between the 
Q^ and the E^ scales in more detail, the 
following cut scenario is adopted. All 
jets lying well within the detector accep- 
tance — 1.5 < rjiab < 2.5 are ordered in 
energy. Then the two highest Et jets 
are required to fulfil: Eti > -E't2 > 
5 GeV and En + Et2 > 17 GeV . 
Jets are defined by the inclusive Kt al- 
gorithm |T3|, ^ in the Breit system. 
The kinematic phase space is defined by 
0.2 < y < 0.6 and 5 < Q^ < 600 GeVl 
These cuts are very close to the ones 
used by the HI collaboration in a recent 



extraction of the gluon density [^ 

In a NLO calculation similar results 
are obtained with the inclusive Kt algo- 
rithm and with the cone algorithm with 
a cone size Rcone = 0.7. This is illus- 
trated in Fig. H where the total dijet 
cross section dPadijet/dQ'^dET in two Q^ 
and mean Et = (-Et,i + -St,2)/2 binsQ 
is shown as a function of the cone size 
Rcone used in the cone algorithm^. In 
case of the cone algorithm the scale de- 
pendence is only slightly infiuenced by 
the cone size. Only at large Et small cone sizes lead to improved scale dependencies. For 



Figure 5: Dependence of the double differential di- 
jet cross section d'^adijet/dQ^dET defined by the cone 
algorithm with a cone size of Rcone as calculated by 
NLO QCD (Et scale) in different bins of Q"^ and of 
the mean Et in the Breit frame. The shaded band in- 
dicates the variation of the NLO prediction for A- li^^^^ 
and 1/4 -/x^Q^. Overlay ed as hatched band is the cross 
section defined by the inclusive Kt algorithm always 
for R = 1. Jets are defined in the Breit frame with 
Eti > Et2 > 5 GeV, Eti + ^T2 > 17 GeV and 
-1 < Vlab < 2.5. 



^Here and in the following figures, d? adijet / {dQ'^ dEr) is the total dijet cross section integrated over a Q^ 
and Et bin and divided by the bin width. 

^For the Kt algorithm the distance parameter R is always set to 1. For an exact definition see ref. pq|. 
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Figure 6: Dijet cross section as a function of Q"^ as calculated by NLO (line) and LO (dashed 
line) QCD using Q^ (a) or the mean E^ (c) as renormalisation scale fi^^.^- ^^^ shaded band 
indicates the variation of the NLO prediction for 4 ■ fj^nom '^"'^ 1/4 ■ f^nom- The ratio of the NLO 
to the LO calculation is shown in (b) and (d). Jets are defined by the inclusive Kt algorithm 
in the Breit frame with Eti > -Et2 > 5 GeV and —1 < rjiab < 2.5. Shown are cross sections 
for Eti + Et2 > 17 GeV and En + Et2 > 40 GeV. 



Rcone = 0.7 both jet algorithms exhibit the same scale dependence. However, the inclusive Kj- 
algorithm is defined in a less ambiguous way and is better suited for NLO comparisons with 
data |]I9|, EOl |21|. Therefore the inclusive K^ jet algorithm is the preferred choice and is used 



in the following. 

In Fig. ^ [Q^ scale) and §c {E^ scale) the inclusive dijet cross section is shown as a function 
of Q2 for Eti + Et2 > 17 GeV and Eti + Et2 > 40 GeV. At low Q^ ^here E^ is higher than 
Q^ the corresponding cross section calculated with E'^ as scale is lower, since as is probed at a 
higher scale. Starting at Q^ > 40 GeV^ the dijet cross section is of comparable size for both 
scales. For both Et cuts, the scale dependence is largest at lowest Q^. The i^-factor shows a 
similar behaviour (see Fig. ^b and Fig. ^d). When using Q^ as scale, at Q^ ~ 10 GeV^ for 
Eti + Et2 > 17 GeV (^ri + Et2 > 40 GeV) the LO cross section is about a factor of 2 (1.75) 



below the NLO cross section. Only at the largest Q^ of 600 GeV^ the NLO correction is small, 
i.e. about 15% (5%). When using E"^ as scale (see Fig. ||c and Fig. |^d), the scale dependence 
is generally reduced and the K-factor is smaller. 
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Figure 7: Double differential dijet cross section d'^adijet/ {dQ'^dEr) as a function of ^ren for 
different bins of Q^ and the mean E\. Shown are NLO and LO QCD predictions for Q^ and 
the mean E^ as renormalisation scale /x^g„ = C,ren ■ fJ'nom- J^^^ ^'^s defined by the inclusive 
Kt algorithm in the Breit frame with Eti > Et2 > 5 GeV, En + Et2 > 17 GeV and 
-1 < Vlah < 2.5. 



The explicit dependence of the double differential dijet cross section on the scale factor ^ren 
is illustrated in Fig. ^ for three different bins in Q^ and mean E^- Shown are the NLO and 
LO cross sections for Q^ and E^ as renormalisation scale. 

For low Q2 (5 < g2 < 10 GeV^) and low Et (8.5 < Et < 13 GeV) the NLO cross section 
exhibits a strong dependence on ^ren as is shown in Fig. |^a. It varies by a factor 2.5 for 
1/4 < C,ren < 4 whcu Q^ IS used a renormalisation scale. For the Et scale this dependence is 
weaker (factor 1.6 for 1/4 < ^ren < 4), but still large. It is interesting that the LO cross section 
is more stable. It falls only by a factor of 1.5 for both the Q^ and the E^ scale. When the Q^ 
scale is replaced by the E^ scale, the cross section diminishes by a factor of 1.6 (1.2) in case 



of NLO (LO). The K-factors are large, i.e. about a factor of 2 — 2.5, for both scales. The E^ 
scale leads to a slightly smaller K-factor. 

Keeping Et fixed at 8.5 < Et < 13 GeV and increasing Q^ to 200 < Q^ < 300 GeV^ 
gives the configuration displayed in Fig. |^. Both choices of the renormalisation scale result in 
similar NLO and LO cross sections. The cross section is now slightly bigger in case of the E^ 
scale. The NLO cross section still reveals a distinct dependence on C,ren, but this dependence 
is now weaker than in the case of the LO cross section. The K-factors are reduced compared 
to the low Q^ result, but are still large (factor 1.5). At no value of C,ren the NLO and the LO 
cross section are similar. 

Fig. 0d presents the results for 5 < Q^ < 10 GeV^ and 30 < Et < 60 GeV. The dependence 
of the NLO cross section on C,ren is much reduced for both choices of scales. The LO cross section 
exhibits a much stronger dependence. For large Et and 5 < Q^ < 300 GeV^ the LO cross 
section is higher than the NLO cross section when Q"^ is used as scale. In case of the Q^ scale 
around C,ren ~ 2 the K-factor is 1. This point moves to higher C,ren values as Q^ increases. In 
case of the E^ scale such a point is at very low ^ren values at low Q^ and is around ^ren = 1/2 
for 35 < g2 < 300 GeV^ 

If both Q^ and Et are large, the NLO cross sections are stable for both scale choices (see 
Fig. 1^). Moreover, the magnitude of the cross section is very similar for both scales. The LO 
cross section is steeply falling as C,ren increases. The K-factor is 1 for ^ren = 1/2 {C,ren = 4) in 
case of the E^ (Q^) scale. 

It is interesting to note that the variation of the NLO cross section with the scale factor ^^en 
as well as the K-factor strongly depend on the position of the jets with respect to the proton. 
This is illustrated in Fig. ^ where the dijet cross section as a function of pseudo-rapidity of 
the forward jet in the laboratory frame rjfyjd^iab is shown for three Q^ bins. Here E"^ is used 
as renormalisation scale. In configurations where both jets are backward, i.e. small values 
of Vfwd,iab^ small scale dependencies and K-factors around unity are found. However, as the 
forward jet moves towards the proton direction increasingly large scale factors are found. For 
5 < g2 < 10 GeV^ the NLO correction is 600% and even for 200 < Q^ < 300 GeV^ the 
NLO correction is still 300% ! Also the scale dependence of the NLO cross section gradually 
increases from about 30% for small r]f^d,iab to about 70% for large rifyjd,iab ^t low Q^. Although 
somewhat smaller, the same effect is visible at large Q^. These large K-factors suggest that 
NNLO corrections will be important. In the specific context of forward jet production studied 
in a more restricted phase space, this has already been demonstrated p2| . 



5 Conclusions 

The behaviour of fixed order calculations in dijet production at HERA has been investigated. 
A cut scenario where the transverse energy of both jets are required to be above the same 
threshold has to be avoided to get agreement of the NLO calculation with the data. Dijet rates 
for 5 < Q^ < 100 GeV^ can be described by a NLO calculation when using the E^ or the Q^ 
scale. The E^ scale leads to smaller scale dependencies. 

The interplay between the virtuality of the photon Q^ and the mean transverse energy of 
the dijet system in the Breit frame E^ has been studied in detail using dijet cross sections for 
5 < g2 < 600 GeV^ and 8.5 < Et < 60 GeV as observable. 
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Figure 8: Dijet cross section as a function of the pseudo-rapidity of the forward jet r]fwd,iab fof 
different Q^ bins. Shown are NLO and LO QCD predictions for the mean E^ as renormalisation 
scale ^nom- The shaded band indicates the variation of the NLO prediction for 4 ■ fin^m '^'^'^ 
l/4-/i^Q^. Jets are defined by the inclusive Kt algorithm in the Breit frame with Eti > Et2 > 5 
GeV, Eti + Et2 > 17 GeV and -1 < r]iab < 2.5. 

Only if both Q^ and E^ are large, the pQCD seems to make reliable predictions. The NLO 
calculation is independent of the renormalisation scale factor ^ren while the LO calculation 
strongly depends on it. Moreover, a point can be found where the NLO and the LO cross 
section have approximately the same size. If either Q^ or E^ is large and the other scale small, 
the NLO is more stable than the LO calculation, but exhibits nevertheless a dependence on 
^ren- This is more pronounced, if Q^ is large and E^ small than in the opposite case where E^ is 
large and Q^ small. When both scales are relatively small, large scale dependencies are found. 
In this case a LO calculation seems to be more stable than the NLO calculation. Moreover, 
large NLO corrections are found. 

Generally, the use of E^ as renormalisation scale is favoured over Q^ , since scale dependen- 
cies are less pronounced and NLO corrections are smaller. 
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